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REDUCTION  OF  FREE-EDGE  STRESS  CONCENTRATION 

P.  R.  Heyliger  and  J.  N.  Reddy 
Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 

Abstract.  The  paper  describes  a  quasi-three  dimensional 
formulation  and  associated  finite  element  model  for  the  stress  analysis 
of  a  symmetric  laminate  with  free-edge  reinforcement.  Numerical  results 
are  presented  to  show  the  effect  of  the  reinforcement  on  the  reduction 
of  free-edge  stresses.  It  is  observed  that  the  interlaminar  normal 
stresses  are  reduced  considerably  more  than  the  interlaminar  shear 
stresses  due  to  the  free-edge  reinforcement. 

1.  Introduction 

Free  edge  effects  in  laminated  composites  have  been  studied  by 
numerous  investigators  [1-14].  These  studies  indicate  the  existence  of 
large  interlaminar  stresses  near  the  free  edge  of  composite  laminates 
under  inplane  loads. 

A  proposed  method  of  reducing  the  intensity  of  these  large  stresses 
is  to  reinforce  the  stiffness  at  the  free  edge  by  wrapping  a  lamina  or 
collection  of  laminae  around  the  free  edge,  forming  a  cap.  The  cap  is 
intended  to  reduce  the  large  displacement  gradients  near  the  free  edge, 
thereby  decreasing  the  interlaminar  stresses  and  reducing  the  proba¬ 
bility  of  delamination.  Experimental  studies  indicate  that  a  rein¬ 
forcing  cap  is  a  viable  method  of  maintaining  laminate  strength  [15]. 

The  purpose  of  this  study  is  to  determine  the  effect  of  a  rein¬ 
forcing  cap  on  the  stress  distribution  in  a  symmetric  composite  laminate 
under  uniform  axial  extension.  Parametric  studies  are  performed  by 
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varying  the  thickness,  length,  and  lamination  scheme  of  the  cap. 
Numerical  results  are  presented  for  several  representative  test  cases. 


2.  Problem  Formulation 

The  initial  work  of  Pipes  and  Pagano  [1]  used  the  three-dimensional 
elasticity  equations  to  develop  an  assumed  displacement  field  for  a 
symmetric  laminate  under  uniform  axial  extension  (Figure  1).  Using  the 
stress-strain  relations  for  a  monoclinic  material  and  the  strain-dis¬ 
placement  relations  for  a  linearly  elastic  material,  symmetry  and  anti¬ 
symmetry  conditions  may  be  imposed  to  give  the  following  displacement 
field: 


u  =  kx  +  U(y,z) 

V  =  V(y,z)  (1) 

w  =  W{y,z) 

Since  all  stresses  are  independent  of  the  axial  coordinate  x,  the 
equilibrium  equations  reduce  to: 


T  +  T  =  0 

xy,y  xz,z 

O  +  T  =0 

y.y  yz,2 

T  +0  =0 

yz,z  z,z 


(2) 


Substitution  of  the  displacement  field  and  stress-strain  relations  into 
the  equilibrium  equations  yields  the  following  governing  equations: 


^66^,yy  ^55^, zz  W,yy  ^45'',zz  ^45^^,yz  "  ° 

^26^yy  "  C45^zz  "  ^22\y,  "  V,zz  (^23  "  C44)W.yz  =  «  <3) 


^^45  46^^,yz  ^^44  ^23^'',yz  ^44'^,yy  ^33*^, zz  "  ° 


These  coupled,  second  order  partial  differential  equations  represent  a 
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u  _(y,o)  =  0 
V  Jy,o)  =  0 


three  variable  problem  in  a  two-dimensional  domain,  and  are  not  solvable 
by  closed  form  methods. 

Applying  symmetry  conditions  about  the  x-z  and  y-z  planes  reduces 
the  problem  domain  to  one  quadrant  of  the  laminate  and  gives  the  fol¬ 
lowing  boundary  conditions: 

U(0,z)  =  0 

V(0,z)  =  0  (4) 

W(y.O)  =  0  ° 

In  addition  to  the  displacement  boundary  conditons,  the  stresses 

may  be  specified  on  the  outer  surfaces  of  the  laminate.  Along  the  top 

surface,  where  the  surface  normal  is  parallel  to  the  z  axis, 

g  =  X  =  T  =0  (5) 

z  yz  xz 

and  along  the  outer  edge,  where  the  surface  normal  is  parallel  to  the  y 
axis, 

a  =  T  =  T  =0  (6) 

y  xy  yz 

With  the  problem  completely  posed,  it  appears  that  the  solution  is 
trivial  (i.e.,  all  of  the  displacements  are  zero  at  every  point  in  the 
domain)  since  the  differential  equations  and  boundary  conditions  are  all 
homogeneous.  This  is  naturally  not  the  case,  but  the  source  of  the 
forcing  function  is  not  obvious  and  will  be  shown  in  sequel. 

The  formulation  summarized  above  may  easily  be  extended  to  include 
the  effects  of  a  cap.  Indeed,  the  governing  equations  and  boundary 
conditions  given  above  will  not  change  for  a  capped  laminate,  with  the 
exception  of  the  change  in  material  properties  in  the  cap.  The  natural 
boundary  conditions  given  in  Eq.  (6)  still  apply  to  the  small  surface 


area  introduced  by  the  cap  whose  normal  is  in  the  negative  y 
direction.  The  inclusion  of  the  cap  material  must  not  violate  any  of 
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the  assumptions  used  to  develop  the  governing  equations  and  boundary 
conditions  (see  [1]).  For  this  reason,  only  cross-ply  lamination 
schemes  may  be  used  for  the  cap  since  the  stress-strain  relations  for  a 
monocllnic  material  require  symmetry  about  the  midplane  of  the  lami¬ 
nate.  The  domain  of  the  new  problem  Is  shown  In  Figure  2  with  the 
necessary  geometric  variables  labeled. 


3.  Method  of  Solution 
3.1  Variational  Formulation 

The  finite  element  method  Is  used  to  obtain  an  approximate  solution 
to  the  three  differential  equations  and  associated  boundary  conditions 
specified  above.  Using  the  Galerkin  type  formulation,  the  three  equa¬ 
tions  are  weighted  with  the  variations  6U,  sV  and  6W  of  the  displace¬ 
ments  U,  V,  and  W,  respectively,  and  are  Integrated  over  the  problem 
subdomain  (element)  R®  and  set  equal  to  zero  (see  [16]).  Considering 
the  first  equation,  this  gives 


J  .  (C35  .  C,5)W_^JdR  =  0 


(7) 


Applying  Green's  theorem  gives 


°  ^  -  C26^y^^y  ^  W,z^^z 

"  '45W,y^^z  "  C36^z^^y)dR 
^  ^55^,z"z  ^26'',y"y  W,z"z 


^45'^,y"z  ^ 


(8) 


In  order  to  Investigate  the  entries  In  the  surface  Integral,  we 
first  note  that  the  stress-strain  relations  for  a  monocllnic  material 
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give 


and 


xy 


(9) 

le'*  *'26'',y  *  *  *'66'^, y' 

(10) 

Comparing  the  entries  of  these  equations  with  those  given  in  Eq.  (8),  we 
may  immediately  rewrite  the  surface  integral  as 


s® 

Once  the  element  assembly  has  been  completed,  the  stresses 

T  and  T  will  cancel  internally  everywhere  in  the  domain  except  on 
xy  xz 

the  outer  surfaces,  where  they  are  specified  to  be  zero  (t  along  the 

xy 

outer  edge  and  t  along  the  upper  surface).  Since  the  material 
properties  may  be  different  for  each  element,  the  term  in  the  surface 
integral  involving  the  initial  strain  k  will  be,  in  general,  nonzero 
between  layers  and  at  the  outer  surfaces  of  the  laminate.  Hence  the 
surface  integral  is  the  source  of  the  forcing  function  and  will  result 
in  a  non-trivial  solution  to  the  governing  equations.  The  other  two 
equations  are  treated  in  similar  fashion  and  yield  similar  results, 
although  the  stresses  in  the  surface  integral  and  the  material 
properties  are  different  for  each  equation. 


3.2  Finite  Element  Equations 

To  implement  the  finite  element  method,  we  first  approximate  the 
displacements  and  their  variations  over  a  given  element  by  assuming  them 
to  be  of  the  form  (see  [16]) 

n 

U  =  z  N.U.  6U  =  N. 

J  J  1 
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n 

V  =  i  N.V  6V  =  N.  (12) 

j=l  J  J  1 


n 

W  =  r  N.W,  6W  =  N. 

j=l  J  J  ^ 

where  Nj  are  the  interpolation  functions  and  n  is  the  number  of  nodes  in 
the  element.  Substituting  these  approximations  into  the  variational 
form  of  the  governing  equations  results  in  the  finite  element  equations 
given  below  for  a  typical  element: 


Ik'^I 

{u}' 

’{F^}  ' 

IK^^l 

[K^^l 

{V} 

"  =  ^ 

{F^}  r 

31 1 

[K^^l 

[K^^j 

{w} 

(13) 


Assembling  the  element  equations  for  the  entire  domain  results  in  the 
simultaneous  equations 

[Kl{a}  =  {F}.  (14) 
The  forms  of  the  local  stiffness  matrix  and  force  vectors  are  given  in 
the  appendix. 


4.  Numerical  Results 

4.1  Interlaminar  Stresses  in  Uncapped  Laminates 

Before  investigating  the  effects  of  a  cap  on  a  symmetric  composite 

laminate,  analyses  were  performed  on  several  uncapped  laminates  to 

ascertain  the  accuracy  of  the  formulation  developed  earlier  in  this 

paper  by  comparing  the  values  of  the  interlaminar  stresses  a  and  x 

z  xz 

with  those  reported  by  other  investigators  for  similar  laminates.  The 
4-ply  laminates  [90“/0°]g  and  (±45° 1^  have  been  studied  extensively  (see 
[9]  and  [14])  in  the  analysis  of  the  free-edge  stress  problem  and  were 
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studied  here  for  purposes  of  comparison.  A  width  to  thickness  ratio  of 
10  was  used  for  these  example  laminates,  and  both  were  subjected  to  an 
initial  strain  of  0.001  in. /in.  The  following  material  properties  were 
usd: 

=  20.0  X  10®  psi 
E2  =  2.1  X  10®  psi 
E^  =  2.1  X  10®  psi 

®12  "  ^ 

G23  =  0.85  X  10®  psi 


^12  ~ 

vi3  =  0.21 

'’23  " 

The  finite  element  mesh  used  to  model  the  upper  right  hand  quadrant 
of  the  two  laminates  contained  160  linear  quadrilateral  elements  and  189 
nodes.  All  stresses  were  computed  at  the  centers  of  the  elements  of  the 
top  ply  that  were  closest  to  the  interface. 

Figures  3-a,b  show  the  interlaminar  stresses  and  plotted 
across  the  width  of  the  [±45° Ig  laminate.  Figure  3-c  shows  the 
distribution  of  the  interlaminar  stress  across  the  width  of  the 

[90°/0°]5  laminate. 

The  results  for  the  [±45° laminate  compare  favorably  with  values 
reported  graphically  in  earlier  works  (e.g..  Ref.  [14]).  Numerical 
values  of  the  stresses  in  both  laminates  are  in  excellent  agreement  with 
the  results  developed  in  [9]  using  an  identical  mesh. 
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4.2  Interlaminar  Stresses  in  Capped  Laminates-Prel iminary  Remarks 

In  order  to  determine  the  effect  of  cap  thickness,  length,  and 
lamination  scheme  on  the  reduction  of  the  interlaminar  stresses,  several 
different  cap  configurations  were  added  to  an  8-ply  symmetric  laminate 
(+30“/-30°/90°/90°)5.  thickness  of  each  ply  is  0.005  in.  This 
laminate,  when  strained,  experiences  both  normal  and  shear  interlaminar 
stresses.  The  material  properties  used  were  that  of  T300/5208  graphite- 
epoxy  and  are  given  below. 

=  19.2  X  10®  psi 
=  1.56  X  10®  psi 
E^  =  1.56  X  10®  psi 
=  0.82  X  10®  psi 
=  0*82  X  10®  psi 

^23  ~  ^ 

=  0.24 

^13  =  0.24 
^23  ^ 

The  dimensions  of  the  problem  domain  are,  for  the  upper  right  hand 
quadrant  of  the  laminate,  b  =  1.0  in.  by  h  =  0.02  in.  (see  Fig.  2); 
hence,  the  width  to  thickness  ratio  of  the  laminate  is  50.  The  laminate 
was  subjected  to  an  initial  strain  of  0.001  in. /in.  in  all  analyses. 

The  same  laminate  was  used  in  all  analyses  mentioned  below. 

The  three  different  cap  configurations  used  in  the  study  are 
identical  except  for  the  length  of  the  overlap  on  the  upper  surface  of 
the  laminate.  This  parameter  is  labeled  L  in  Figure  2.  The  three 
different  values  used  were  0.4  in.  (cap  one),  0.15  in.  (cap  two),  and 
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6.25  X  10"^  (cap  three).  Perfect  bonding  is  assumed  to  exist  between 
ail  elements. 

The  lamination  scheme  of  the  cap  must  be  limited  to,  at  the  most,  a 
cross-ply  due  to  the  symmetry  requirements.  Each  cap  configuration  was 
analyzed  using  two  different  lamination  schemes;  one  with  all  the  cap 
fibers  perpendicular  to  the  x-axis,  and  the  other  with  the  cap  fibers 
parallel  to  the  x-axis.  The  fibers  of  the  first  cap  are  assumed  to  bend 
at  right  angles  at  the  upper  right  hand  corner  of  the  cap,  with  the 
fibers  along  the  top  and  edge  parallel  to  the  y-  and  z-axis,  respec¬ 
tively. 

In  addition  to  varying  the  overlap  length  and  lamination  scheme  of 
the  cap,  the  effect  of  cap  thickness  (labeled  T  in  Figure  2)  was  also 
investigated.  The  three  values  used  were  0.02  in.,  0.01  in.,  and  0.005 
in.  The  same  mesh  was  used  for  all  analyses  except  for  the  changes  in 
the  cap  element  dimensions  due  to  the  change  in  cap  thickness.  The 
meshes  used  to  model  the  uncapped  and  capped  laminates  are  shown  in 
Figs.  4a  and  4b,  respectively. 

Before  varying  the  parameters  of  the  cap,  it  is  interesting  to 
observe  the  changes  in  the  interlaminar  stresses  along  the  width  of  the 
laminate  due  to  the  addition  of  a  cap.  The  8-ply  laminate  was  first 
analyzed  using  the  meshes  shown  in  Figs.  4a  and  4b.  Both  meshes  were 
then  slightly  refined  near  the  free  edge  and  near  the  interface  and 
midplane  elements  obtain  a  better  approximation  of  the  interlaminar 
stress  distributions  for  the  uncapped  and  capped  laminates.  An  exploded 
view  of  the  capped  portion  of  the  refined  mesh  is  shown  in  Fig.  4c.  The 
refined  mesh  for  the  uncapped  laminate  is  not  shown. 
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Figures  5a-b  show  the  stress  distributions  of  o  and  t  plotted 
against  the  y  coordinate  of  the  laminate  for  a  capped  (L  =  0.15  in.,  T  = 
0.005  in.,  fibers  perpendicular  to  x-axis)  and  uncapped  laminate  for  the 
two  different  meshes.  The  values  for  were  taken  from  elements  close 
to  the  laminate  mid-plane,  and  the  values  for  were  taken  near  the 
interface  of  the  +30°  and  -30°  plies.  Only  the  stresses  in  the  region 
near  the  free  edge  are  plotted  since  the  interlaminar  stresses  dissipate 
at  interior  regions  of  the  laminate.  The  solid  lines  represent  the 
stresses  in  the  uncapped  laminate,  and  the  dashed  lines  represent  the 
stresses  in  the  capped  laminate.  All  curves  shown  in  this  report  were 
plotted  using  stresses  taken  from  the  center  of  the  elements  and  were 
smoothed  by  a  spline  routine. 

These  four  figures  show  that  a  capped  laminate  exhibits  reduced 
interlaminar  stresses  near  the  free  edge  when  compared  with  the  uncapped 
laminate.  The  interlaminar  stress  is  reduced  by  a  sizeable  amount  at 
the  data  point  closest  to  the  free  edge,  while  the  reduction  of  the 
stress  T  is  smaller  but  still  noticable.  The  x  stress  distribution 

aZ  X2 

is  graphically  similar  for  both  the  capped  and  uncapped  laminates,  with 
the  stresses  in  the  capped  laminate  slightly  lower  than  the 
corresponding  stress  at  the  same  location  in  the  uncapped  laminate. 

The  distribution  is  significantly  different  for  the  capped  and 
uncapped  laminates,  showing  an  oscillatory  type  of  behavior  near  the 
free  edge  of  the  capped  laminate  before  converging  to  zero  further 
inside  the  laminate.  The  length  of  the  edge  effect  is  greatly  increased 
by  the  presence  of  a  cap. 

Figures  5a-d  show  that  the  interlaminar  stresses  in  the  capped 
laminates  do  not  approach  a  finite  value  near  the  free  edge.  It  also 
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appears  from  these  figures  that  the  slopes  of  the  curves  appear  to  be 
Increasing  near  the  edge  of  the  laminate.  Hence  It  Is  unlikely  that  the 
stress  singularity  has  been  eliminated.  However,  the  figures  do  show 
that  the  magnitudes  of  the  free  edge  stresses  are  reduced  somewhat  by 
the  inclusion  of  a  cap  on  a  symmetric  laminate  under  uniform  axial 
extension,  and  indicate  that  a  reinforcing  cap  may  reduce  the 
probability  of  delamination  in  such  a  laminate. 

Subsequent  analyses  will  be  concerned  with  the  effect  of  cap 
parameter  variation,  described  earlier  in  this  section,  n  interlaminar 
stress  reduction  near  the  free  edge. 

4.3  Results  and  Discussion 

The  results  of  the  analyses  are  summarized  in  Figs.  6  8.  The 

stresses  in  the  uncapped  laminate  are  represented  by  a  solid  line  in  all 
of  the  graphs.  The  maximum  normal  stress  occurs  near  the  midplane  of 
the  laminate,  and  the  maximum  shear  stress  occurs  near  the  +30V-30° 
interface.  Since  the  stresses  are  calculated  at  the  element  centers, 
the  plots  do  not  show  the  stresses  exactly  at  the  midplane  or  the  lami¬ 
nate  edge. 

Figures  6-8  show  that  the  stresses  and  for  the  laminates 
with  caps  one  and  two  are  nearly  identical  for  all  parameter  varia¬ 
tions.  The  laminate  with  cap  three,  however,  experiences  larger  stress 
reductions  when  compared  with  the  other  capped  laminates.  Thus  it 
appears  that  the  overlap  length  of  the  cap,  once  past  a  certain  point, 
has  little  if  any  effect  on  the  performance  of  the  cap. 

Although  both  of  the  lamination  schemes  studied  give  nearly  iden¬ 
tical  distributions  for  the  shear  stress  the  cap  with  the  fibers 
parallel  to  the  x-axis  induces  greater  reductions  in  the  normal 
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stress  a^.  The  thinner  caps,  in  general,  create  larger  normal  stress 
reductions  for  the  caps  with  perpendicular  fibers,  while  the  caps  with 
the  fibers  parallel  to  the  x-axis  show  more  variability. 

One  important  aspect  of  the  feasibility  of  a  cap  is  the  stress 
arising  in  the  cap  itself.  The  purpose  of  the  reinforcing  cap  is  to 
provide  stiffness  through  the  thickness  of  the  laminate  at  the  free 
edge,  where  the  weakness  of  the  matrix  material  in  tension  could  result 
in  delamination.  Capping  the  laminate  may  indeed  reduce  the  inter¬ 
laminar  stresses  in  the  laminate  itself,  but  could  cause  large  stresses 
within  the  cap.  Referring  to  Fig.  5a  and  noting  that  the  asterisks 
represent  the  stress  values  in  the  cap  elements,  we  may  observe  the 
widely  varying  values  of  occurring  in  the  cap.  This  may  indicate  one 

disadvantage  of  using  a  cap  comprised  strictly  of  fibers  parallel  to  the 
x-axis,  since  this  type  of  cap  would  be  just  as  weak  as  the  laminate  in 
the  z  direction.  A  cap  with  a  cross-ply  lamination  shceme  could  prove 
beneficial  to  provide  stiffness  and  strength  in  both  x  and  z  direc¬ 
tions.  Typical  stress  distributions  for  two  different  cross-ply  schemes 
are  shown  in  Fig.  9.  This  figure  represents  a  two-ply  cap  with  the 
dotted  lines  and  dashed  lines  representing  parallel  fibers  on  the  outer 
ply  and  inner  ply  of  the  cap,  respectively. 

5.  Conclusions 

The  analyses  completed  in  this  study  show  that  a  reinforcing  cap 
can  reduce  the  interlaminar  stresses  near  the  free  edge  in  a  symmetric 
composite  laminate  under  uniform  axial  extension.  The  overlap  length  of 
the  cap  makes  little  difference  in  the  effectiveness  of  most  caps  that 
may  realistically  be  fabricated.  Lamination  schemes  are  relatively 
equivalent  in  reducing  shearing  stresses,  but  may  be  of  greater  concern 
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regarding  normal  stress  reduction.  Caps  of  all  thicknesses  reduced  the 
stresses  in  the  laminates,  but  no  conclusions  concerning  an  ideal  thick¬ 
ness  may  be  readily  drawn  from  the  results  reported  herein. 

Although  a  single  laminate  was  used  for  purposes  of  comparison  in 
this  study,  the  general  trends  are  the  same  for  other  laminates  with 
large  interlaminar  stresses.  It  is  unlikely  that  one  cap  configuration 
is  ideal  for  all  laminates,  and  specific  caps  must  be  designed  to 
account  for  the  characteristics  and  behavior  of  particular  laminates. 
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APPENDIX 

Coefficients  of  the  stiffness  matrix  in  Eq.  (13) 
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Figure  1  Coordinate  system  for  the  laminate  under  consideration 
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Figure  2  Domain  of  the  problem  (upper  right  quadrant  of 
a  symmetric  capped  laminate) 
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Figure  3  (b)  x  distribution  across  width  of  laminate 
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Figure  3  (c)  distribution  across  width  of  laminate 


(a)  Original  mesh  for  uncapped  laminate 


(b)  Original  mesh  for  capped  laminate 


(c)  Edge  section  of  refined  mesh  for  capped  laminate 


Figure  4  Finite  element  meshes  used  to  model  the  uncapped  and  capped  laminates 
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Figure  5  (a)  distribution  for  capped  and  uncapped  laminate  (mesh  4a-b) 
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Figure  5  (b)  t  distribution  for  capped  and  uncapped  laminates  (mesh  4a-b) 
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Figure  5  (d)  t  distribution  for  capped  and  uncapped  laminate  (refined  mesh) 
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Figure  6(a) 


az  distribution  along  laminate  edge  -  cap  1 
(  fibers  parallel  to  the  x-axis  ) 
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Figure  6  (b)  distribution  along  laminate  edge  -  cap  2 
(  fibers  parallel  to  the  x-axis  ) 
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Figure  7(a)  distribution  along  laminate  edge  -  cap  1 
(  fibers  perpendicular  to  the  x-axis) 
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Figure  7(b)  distribution  along  laminate  edge  -  cap  2 
(  fibers  perpendicular  to  the  x-axis  ) 
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Figure  7(c)  distribution  along  laminate  edge  -  cap  3 
(  fibers  perpendicular  to  the  x-axis  ) 
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T  distribution  along  laminate  edge  for  cap  3 
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Figure  8 (a) 
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T  distribution  along  laminate  edge  for  cap  3 
xz 

(  fibers  perpendicular  to  the  x-axis  ) 
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Figure  9(a)  distribution  along  laminate  edge  for  cross-ply  cap  2 
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Figure  9(b) 


T  distribution  along  laminate  edge  for  cross-ply  cap  2 
xz 


UNCLASSIFIED 


security  cl  ASSI  FtCATlON  OF  THIS  PAGE  (When  Duta^F.ntf^rtui) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT'NUMBER 

VPl-E-85-2 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

«.  TITLE  C«nd  Subr/(/«) 

REDUCTION  OF  FREE-EDGE  STRESS  CONCENTRATION 

5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

Interim 

6.  performing  org.  report  number 
Tech.  Report  No.  ONRM-85-1 

7.  AUTHORf*; 

P.  R.  Heyliger  and  J.  N.  Reddy 

B.  CONTRACT  OR  GRANT  NUMBERf.J 

N00014-84-K-0552 

9.  PERFORMING  ORGANI  Z  ATION  N  AM  PI  AND  ADDRESS 

Virginia  Polytechnic  Institute  and  Slate  Universi 
Blacksburg,  Virginia  240b  1. 

10.  PROGRAM  element,  PROJECT,  TASK 
AREA  «  WORK  UNIT  NUMBERS 

ty 

NR-064-727/5-4-84  (430) 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

Mechanics  Division  (Code  430) 

800  N.  Quincy  St.,  Arlington,  VA  22217 

12.  REPORT  DATE 

January  1985 

13.  NUMBER  OF  PAGES 

35 

14  monitoring  agency  name  ft  ADDRESSf//  <f//ferenf  from  CnntrnfUnfi  Offirn) 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

'iSo.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 

16  DISTRIBUTION  STATTMFNT  (at  thi» 

This  document  has  been  approved  for  public  release  and  sale;  distribution 
uul Lmi ted 

17.  distribution  statement  (of  the>  Abstract  entered  In  Block  70,  If  different  from  Report) 

10.  SUPPLEMENTARY  NOTES 

Thp  report  is  submitted  as  a  paper  for  publication  to  Journal  of  Applied  Mech. 

The  paper  is  also  accepted  as  a  student  paper  for  presentation  at  the  AIAA 

2(ith  Structures,  Structural  Dynamics  and  Materials  Conference,  Orlando,  FL. 

19  KfY  WORDS  ofi  revi^t.ie  nitlr  If  not  etiunry  anti  hlrntlfv  tty  />/<•«  A  nii/ti/mr) 

Composite  laminates,  finite  element  analysis,  free-edge  reinforcement,  free- 
edge  stresses,  interlaminar  stresses,  quasi-three-dimensional  analysis, 
stress  reduction,  symmetric  laminates. 

20  ABSTRACT  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

A  quasi-three-dimensional  formulation  and  associated  finite-element  analysis 
for  the  stress  analysis  of  symmetric  laminates  with  free-edge  reinforcements 
are  presented.  Numerical  results  are  presented  to  show  the  effect  of  the 
reinforcement  on  the  reduction  of  free-edge  stresses.  It  is  observed  that  the 
normal  stresses  are  reduced  considerably  more  than  the  interlaminar  shear 

stresses  due  to  the  free-edge  reinforcement. 

DD  ,  ^°r73  1473  EDITION  OF  1  NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


bLCURMY  CL  AbSIf  K:ATI0N  Of  TmS  F’AGF  (Wi>**n  I>»f» 


VIRGINIA  TECH  CENTER  FOR 
COMPOSITE  MATERIALS  AND  STRUCTURES 


The  Center  for  Composite  Materials  and 
Structures  is  a  coordinating  organization  for 
research  and  educational  activity  at  Virginia 
Tech.  The  Center  was  formed  in  1982  to 
encourage  and  promote  continued  advances  in 
composite  materials  and  composite  structures. 
Those  advances  will  be  made  from  the  base  of 
individual  accomplishments  of  the  thirty-four 
founding  members  who  represent  ten  different 
departments  in  two  colleges. 

The  Center  functions  by  means  of  an 
Administrative  Board  which  is  elected  yearly. 
The  general  purposes  of  the  Center  include: 

•  collection  and  dissemination  of  informa¬ 
tion  about  composites  activities  at  Virginia 
Tech, 

•  contact  point  for  other  organizations  and 
individuals, 

•  mechanism  for  collective  educational  and 
research  pursuits, 

•  forum  and  mechanism  for  internal  inter¬ 
actions  at  Virginia  Tech. 

The  Center  for  Composite  Materials  and 
Structures  is  supported  by  a  vigorous  program 
of  activity  at  Virginia  Tech  that  has  developed 
since  1963.  Research  expenditures  for  investiga¬ 
tions  of  composite  materials  and  structures  total 
well  over  five  million  dollars  wi*  yt  arly 
expenditures  presently  approachi.  illion 

dollars. 


Research  is  conducted  in  a  wide  variety  of 
areas  including  design  and  analysis  of  compo¬ 
site  materials  and  composite  structures, 
chemistry  of  materials  and  surfaces,  characteri¬ 
zation  of  material  properties,  development  of 
new  material  systems,  and  relations  between 
damage  and  response  of  composites.  Extensive 
laboratories  are  available  for  mechanical 
testing,  nondestructive  testing  and  evaluation, 
stress  analysis,  polymer  synthesis  and  character¬ 
ization,  material  surface  characterization, 
component  fabrication  and  other  specialties. 

Educational  activities  include  eight  formal 
courses  offered  at  the  undergraduate  and 
graduate  levels  dealing  with  the  physics, 
chemistry,  mechanics,  and  design  of  composite 
materials  and  structures.  As  of  1982,  some  33 
Doctoral  and  37  Master’s  students  have 
completed  graduate  programs  and  several 
hundred  Bachelor-level  students  have  been 
trained  in  various  aspects  of  composite 
materials  and  structures.  A  significant  number 
of  graduates  are  now  active  in  industry  and 
government. 

Various  Center  faculty  are  internationally 
recognized  for  their  leadership  in  composite 
materials  and  composite  structures  through 
books,  lectures,  workshops,  professional  society 
activities  ar  f  research  papers. 


F  A  ,  V  r  v*,  KS  Ok^  TH  ENTER 


Aerospace  and  Ocean 

Engineering 
Raphael  T.  Haftka 
William  L.  Hallauer,  Jr. 
Eric  R.  Johnson 

Chemical  Engineering 
Donald  G.  Baird 

Chemistry 
James  E.  McGrath 
Thomas  C.  Ward 
James  P,  Wightman 

Civil  Engineering 
Raymond  H.  Plant 

Electrical  Engineering 
loannis  M.  Besieris 
Richard  O.  Claus 


tr*;:  ^  :>cience  and 

ie*  •  . 

•  A. 

c:'  t*-  <1  ^ 

c 

\  -  ’  V  r  V  f 

"  ..  h:.  ■ 

,  nohar 

Alfred  C.  U 
Don  H,  Mor«  i 
Daniel  Post 
J.  N.  Reddy 
Kenneth  L  Reifsnif  ^ 

C.  W.  Smith 
Wayne  W.  Stinchcomb 


Industrial  Engineering 

and  Operations  Research 
Joel  A.  Nachlas 

Materials  Engineering 
David  W.  Dwight 
D.  P.  H.  Hasselman 
Charles  R.  Houska 
M.  R.  Louthan,  Jr. 

Mathematics 
Werner  E.  Kohler 

Mechanical  Engineering 
Norman  S.  Eiss,  Jr. 
Charles  E.  Knight 
S.  W.  Zewari 


inquiries  should  be  directed  to: 

Center  for  Composite  Materials  &  Structures 
College  of  Engineering 
Virginia  Tech 
Blacksburg,  VA  24061 
Phone:  (703)  961-4969 


U21822? 


